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Abstract 

Hypercomplex Fourier transforms are increasingly used in 
signal processing for the analysis of higher-dimensional sig- 
nals such as color images. A main stumbling block for 
further applications, in particular concerning filter design 
in the Fourier domain, is the lack of a proper convolu- 
tion theorem. The present paper develops and studies 
two conceptually new ways to define convolution products 
for such transforms. As a by-product, convolution theo- 
rems are obtained that will enable the development and 
fast implementation of new filters for quaternionic signals 
and systems, as well as for their higher dimensional coun- 
terparts. 
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1 Introduction 

Recently, there has been an increased interest in applying 
hypercomplex Fourier transforms (FTs) in various aspects 
of signal processing where higher dimensional or vector 
signals are used, such as color image processing [29 j [37 ] . 
flow visualization [25l [26] , and even spoken word recogni- 
tion [3]. The main idea behind these applications is the 
representation of a signal (say, a color image) as a pure 
quaternion or as an element of a suitable Clifford algebra 
(see Section[2]for a precise definition). This representation 
is subsequently analysed using a generalisation of the clas- 
sical Fourier transform to a hypercomplex FT, which takes 
into account the multi-dimensional and multi-component 
nature of the signal under consideration. This stands in 
stark contrast to a component based classical analysis, 
sometimes also called marginal analysis. Successful fur- 
ther developments of the hypercomplex approach include 
the design of a color edge filter |42j , as well as other filters 
[45] , construction of FFT methods to compute hypercom- 
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plex FTs [40], etc. 

The main issue that hinders further development of ap- 
plications (and in particular of filter design not based on 
ad hoc assumptions or ideas) is the lack of a suitable con- 
volution theorem. Indeed, in ^lOj it was shown that hyper- 
complex FTs such as the quaternionic FT do not interact 
nicely with the classical convolution product, but rather 
lead to very complicated expressions in the Fourier do- 
main. This means that, up to now, no filter design was 
possible in the Fourier domain, and hence that no fast im- 
plementations as multiplication operators have been ob- 
tained so far. 

The main aim of the present paper is to tackle that 
problem. We will investigate, on theoretical grounds, the 
different possible convolution products that can be defined 
for a wide class of hypercomplex FTs. It turns out that 
two conceptual ways exist to achieve this, both having a 
left and right version. Before explaining these two def- 
initions, let us recall that in the recent literature three 
different approaches to hypercomplex FTs have been con- 
sidered. We can identify them as follows 

• A: Eigenfunction approach 

• B: Generalized roots of —1 approach 

• C: Characters of Spin group approach 

The first approach is mainly studied in [6j [3 [H [151 [IB 
\T7\ [T8] [T9l [20] [21] , and aims at constructing new hyper- 
complex transforms by prescribing eigenvalues to a suit- 
able basis of a Clifford-algebra valued L2 space. The 
choice of suitable eigenvalues implies that there is a huge 
design freedom in this approach. The transforms of this 
class also have a deep connection with quantum mechanics 
and exhibit a very particular underlying algebraic struc- 
ture, namely that of the Lie superalgebra osp(l|2). For a 
recent review from this point of view, we refer the reader 
to Hg. 

The second approach is mainly advocated in (9] [10] and 
boils down to replacing the imaginary unit i in the ex- 
ponent of the ordinary Fourier transform by a generalized 
root of —1, belonging to a Clifford algebra (see e.g. |5^IM] 
for a detailed study of such roots). It encompasses several 
of the hypercomplex FTs often used in applications, such 
as the quaternionic Fourier transform [29], the Sommen- 
Biilow transform [TTJ [46] , the Clifford Fourier transform 
(written without hyphen) introduced in j26j and further 
extended in [J [33] , and the cylindrical Fourier transform 
|8]. Again this approach exhibits a huge design freedom, 
as the set of roots of —1 is very big and as the roots can 
in principle be chosen independently for each application. 
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Finally, a third approach is given in [21 [3], and uses 
the notion of character (or group morphism) to generalize 
the ordinary Fourier transform to the setting of the group 
Spin{3), resp. Spin{5) for direct application in grey scale, 
resp. color image processing. Although the conceptual 
ideas of this third approach are very different from the sec- 
ond approach B, the resulting transforms can nevertheless 
be written as special cases of the transforms in B. For that 
reason we will only focus on the first two approaches in 
the sequel. 

The classical FT and the classical convolution product 
will serve as a guide to define the generalized convolution 
products. Recall that the classical convolution product for 
two functions / and g is defined by 

f*g{x) := / f{y)Tyg{x)dy, 

with Tygix) := g{x — y). The classical FT, defined over 
R™, is given by 

JR'" 

for / g Li(R''"), where {x,y) — J^JLi ^jVj the standard 
inner product. It interacts nicely with the convolution 
product. Namely, one has 

^(/*5)-(2^r/'-^(/m5). (1) 

A first possibility to generalize the convolution product is 
hence obtained by taking the inverse FT of the right-hand 
side of ([U as a definition. This idea was first explored by 
Mustard for the fractional Fourier transform, see [38 . 

Another way to generalize the convolution product is 
obtained by introducing a generalization of the (geomet- 
ric) translation operator Ty. This is also the strategy which 
is e.g. used for the Dunkl transform (see t41,i49,). Let us 
first illustrate the concept in some more detail for the or- 
dinary FT. First we compute the FT of the translation 
over z oi a function /: 

-F(r,/)(?y) = (27r)-^^ e-'^'-^'y^ f{x - z)dx 

= e-*<"'^>(27r)-^ f e-'^^'y^fix)dx 

JR'" 

= e-'^^^y^T{f){y). 

This means that, formally, the ordinary translation is re- 
covered via 

rJ{u)=T-'[e~^<^-y^nf){y))- (2) 

Here, the inverse FT acts on the y variable and yields the 
u variable. 

Now let K{x,y) be the integral kernel of a hypercom- 

plex FT and K{x, y) be the kernel of its inverse J^]^^ ■ 
Upon replacing e~'^^'^^ by the kernel of the hypercomplex 



FT under consideration in formula ([2]), one obtains the 
definition of a generalized translation operator: 

r^f{x) = [ K{Cx)Kiy,O^K{f){m- (3) 

JR™ 

Note that this definition immediately reduces to geometric 
translation (i.e. formula ([2])) when K is the kernel of the 
classical FT. 

We can hence summarize the four possible definitions 
for a generalized convolution product, related to a hyper- 
complex FT Tk'- 

• The first type of definition is inspired by convolution 
for the fractional Fourier transform. It is given by 

f*i9{x) ■.= T^^ (^TK{f)TK{g)){x) 

or 

/*2g(x) ■.= Tk\Tk{9)Tk{!)){x). 

These two versions in general do not coincide due to 
the non-commutativity of the Clifford algebra. How- 
ever, one immediately observes that f *i g = g *2 I ■ 
For symmetry reasons it is nevertheless convenient to 
introduce both versions. 

This definition is clearly interesting as it forces a con- 
volution theorem to hold for the hypercomplex FT 
under consideration. As we will see in the sequel, it 
is possible to derive explicit formulas for these new 
convolutions in terms of classical convolutions (see 
e.g. the subsequent Theorem IH]). 

• Using the generalized translation operator on the 
other hand, one can also define 

/*3 5(^):= / [r^f{^)]9{y)dy 

JR" 

or 

/*4ff(x) -.^ J^^f{y)[T^g{x)]dy. 

Both definitions are again different due to the non- 
commutativity of the Clifford algebra. Moreover, in 
this case there is no immediate relation expressing 
/ *3 g in terms of / *4 g or vice versa. 

This definition is interesting for a different reason: it 
allows to obtain inversion theorems for hypercomplex 
FTs in a way similar to Theorem 8.4 in 1211. 

For hypercomplex FTs, so far only definition f *3 g has 
been applied, in only one particular case, see [21]. In the 
present paper, we will perform a detailed study of the 
four different products mentioned above, for transforms 
belonging to approach A and B. In both cases, we will in- 
troduce specific notations to distinguish between the var- 
ious definitions. An overview of these notations and the 
results obtained in the paper is given in Table [TJ Note 
that the methods developed in the paper can equally be 
extended to other hypercomplex FTs that have not been 
designed yet. 
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The paper is organized as fohows. After some prchni- 
inaries on Chfford algebras and analysis, we introduce in 
Section [3] the most general hypercomplex FT that can be 
considered within the approach A. Its integral kernel takes 
the form of an infinite series in terms of Gegenbauer poly- 
nomials and Bessel functions. We discuss several examples 
and calculate the action of the associated transform on a 
basis of the relevant function space. We also construct 
the inverse of this general transform. In Section 13.21 we 
define the four types of convolution for this general trans- 
form, and study their basic properties. Next, in Section 
13.31 we compute the action of the generalized translation 
operator on the special class of radial functions and prove 
that, for this special set of functions, the new translation 
coincides in many cases with ordinary translation. The 
proof of this result is technical, and the reader may wish 
to skip it during a first reading of the paper. In Section |4l 
we turn our attention to the transforms belonging to ap- 
proach B. We give the general definition and compute the 
eigenvalues and eigenfunctions (thus establishing a con- 
nection with the techniques used in approach A). Then 
we revisit the different convolution definitions in this con- 
text. Both in the case of Mustard convolution (Section 
14. 2p and in the generalized translation operator approach 
(Section 14. 3p . we can now obtain very explicit formulas 
for the new convolutions, which can moreover easily be 
implemented in software packages. In Section 14.41 we re- 
state our results for the special case of the quaternionic 
Fourier transform (qFT) . This is the most well-known hy- 
percomplex FT used in engineering, and hence deserves a 
separate treatment. We end the paper with Table [T] where 
the results are summarized for the two approaches A and 
B. 

2 Preliminaries on Clifford alge- 
bras and analysis 

The Clifford algebra C/o,m over M™ is the algebra gener- 
ated hy ei, i — 1, . . . ,m, under the relations 

CiCj + ejgj =0, i ^ j 

el = -1. 

This algebra has dimension 2"* as a vector space over M. 
It can be decomposed as CZo,m = ffifcLo^^o.m ^i^h CIq „^ 
the space of fc-vectors defined by 

Clo,m ■= span{eii . . . a^^ii < . . . < ik}. 

In the applied literature, Clifford algebras are usually 
called geometric algebras. For a detailed exposition from 
this point of view, including geometric interpretations and 
applications in computer vision, we refer the reader to [24] . 

In the sequel, we will always consider functions / taking 
values in CZo,?n, unless explicitly mentioned. Such func- 
tions can be decomposed as 

m 

C\ . . . Cm/l...m 



with fo, fi, fij, ■ ■ ■ , fi...m all real- or complex- valued func- 
tions on M™. 

The Dirac operator is given by dx ■— X^Jli^^ij^j ^'^^ 
the vector variable by x :— X^JLi ^j^j- The square of the 
Dirac operator equals, up to a minus sign, the Laplace op- 
erator in R™: = — A. For more information regarding 
Clifford analysis, we refer the reader to 122]. 

Denote by V the space of polynomials taking values in 
CZo,m, i.e. 

P ■.^R[xi,...,Xm](E>Clo,m. 

The space of homogeneous polynomials of degree k is then 
denoted by Vk- The space Mk ■= kerdx<^Vk is called the 
space of spherical monogenics of degree k. 

Next we define the inner product and the wedge product 
of two vectors x and y 

m 

{x,y) ■^^ '■>^jyj 

xAy:^ ^^ejCkixjyk - Xkyj). 

]<k 

Remark 1. From now on, we will use the following 
convention for denoting variables, resp. vectors. When 
X = (xi, . . . , Xrn) S R™ is used as a variable we will not 
underline it (in order not to overload notations). When 
it is used as a vector x :— X^JLi ^j^j involving Clifford 
multiplication, we use the underlined notation. 

We introduce two different bases for the space 
iS(]R™) (8)C^o.mj where iS(R™) denotes the Schwartz space. 
Define the functions '4'j,k,e by 

(4) 

V2,+i,M :=i/+'(NPkMfe^-l^/^ 

where j,k G N, {M^^^ e Mk ■ i = I, . . . ,dimMk} is a 
basis for A4k, and L" are the Laguerre polynomials. The 
set {ipj,k,i} forms a basis of S{W^) (8)CZo,m, see [47]. This 
basis is called the Clifford-Hermite basis or the spherical 
basis. Alternatively, define the one-dimensional Hermite 
functions (see e.g. [JS]) by 

= Hfc(x)e-'/2 
for A; e N. Then the set {V'ji j2....jm} with 

and ji, . . . , S N is also a basis of 5(M™) (g>C^o,m, called 
the tensor product or cartesian basis. Both bases interact 
nicely with the ordinary FT. One has 

-^(^,,m) = (-«)-''+'^j,m (5) 

It is interesting to compare this result with the subsequent 
Theorems [T] and [Tj For more information about these two 
bases, we refer the reader to [T2j. 
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3 Convolution products in Ap- 
proach A 

3.1 Definition of tiie transforms, eigen- 
functions, eigenvalues and inverse 

In this section we consider a general kernel of the following 
form 

K{x, y) ^ {A{w, F) + (x A y)B{w, T)) 

^ g^(cota)(|s.|^ + |y|^) 

with 

+00 

fc=0 

+00 
fe=i 

and ak.Pk £ C, z = (|x||y|)/sinQ!, w = {x,y) / {\x\\y\), 
A = (m — 2)/2 and a G [— 7r,7r]. Here, is the Bessel 
function and the Gegenbauer polynomial. We exclude 
the cases where a = or a = ±7r. 

The integral transform associated with this kernel is de- 
fined by 

if) (y) = Po,™ / K{x, y) fix) dx (8) 

with 

Po.,m = (^(1 ~ e-2-))-W2 

and with dx the standard Lebesgue measure on M'". The 
precise form of the kernel in formula ^ is inspired by the 
results obtained in [161 . In particular, it encompasses all 
previously studied kernels in the eigenfunction approach 
A. 

Remark 2. In order to ensure good analytic behavior 
of the transform (0^, we additionally demand that both 
A(w, z) and B(w, z) satisfy polynomial bounds of the fol- 
lowing type 

\Aiw,I)\<c{l + \x\)Hl + \y\y, 
\B{w,^\<c{l + \x\y{l + \y\y, 

with j € N and c a constant. 

This is the case for all subsequent examples, and allows 
us to prove that the transform (0) yields a continuous map 
on S{W") ® Cla.m, by adaptation of the proof in J Iff]. 

3.1.1 Some examples 

In the special case where a — 7r/2, the kernel takes the 
form 

K{x,y) = A{vu,z) + (xAy) B{w,z) (9) 



with 

+00 

A{w, z)=J2 "fc z-^Jk+x{z)C^{w) 

k=0 
+00 

B{w,z) = z-^-^Jk+x{z)C^ll{w) 

k=l 

and z = |a:||y|, w = {x , y) / {\x\\y\) , A = (m - 2)/2. The 
corresponding integral transform is given by 

J^K if) (y) = Pf .m / K{x, y) f{x) dx 

where p^^„, = {2ny^/\ 

Note that the kernel of the classical Fourier transform 
J-", which can equally be expressed as the operator expo- 
nential eT^eT('^-l=^l )^ takes the form ([9]) with 

ak^2^T{\){k + X){-if 

Also the Clifford-Fourier transform (see O H H [21]), a 
generalization of the classical Fourier transform in the 
framework of Clifford analysis, takes this form. It is de- 
fined by the following exponential operator 

with 

r := -^ejekixjdx^: -xtdx^)- 
j<k 

In case of the Clifford- Fourier transform the coeffi- 
cients ak and /S^ in the kernel © take the form: 

afc = 2^-ir(A + l)(i2^+2 + (-!)'=) 

- 2^-ir(A) (fc + A)(^2^+2 - (-1)^-) 
/3fc = -2^^(A + l)(^2^+2 + (-l)'=). 

For the transform 7^+, similar expressions hold. Moreover, 
in [15] , an entire class of kernels of the form © , for partic- 
ular values of the coefficients ak and /3fc, was determined. 
They yield new integral transforms that have the same 
calculus properties (i.e. interaction with the Dirac oper- 
ator) as the original Clifford-Fourier transform, but with 
different spectrum. 

Also for general a, concrete examples have been studied. 
The fractional Fourier transform (see [55]) is a generaliza- 
tion of the classical Fourier transform. It is usually defined 
using the operator expression 

J'a = e^e^*^^~l^l^\ ae[-7r,7r]. 

Recently, a fractional version of the Clifford- Fourier trans- 
form was introduced (see [TB])- It is defined by the follow- 
ing exponential operator 
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The integral kernel of this transform takes the form ([7]) 
with 

2^r(A + 1) 



Proof. This follows from the explicit expression ^ of the 
basis and the identity (see |31j p. 847, formula 7.421, 
number 4 with a = 1]): 



+ 00 



3.1.2 Eigenvalues 



k( if3{k+2X) _ -i^k 



4/3(1-/?) 



Now we calculate the action of the transform ^ on the 
basis (gl) of 5(E'") (g) CZo,m- We start with the following 
auxiliary result, expressing the radial behavior of the in- 
tegral transform. 

Proposition 1. Let G Mk be a spherical monogenic 
of degree k. Let f{x) = /o(|a:|) he a real-valued radial 
function in 5(R™). Further, put ^ — x/\x\ and r/ = y/\y\- 
Then one has, putting /3o = 0, 



TK{fir)Mk){y) 
X 



X + k 

+ 00 



ak — sm a 



2{k + X)' 



□ 

Remark 3. Theorem{l\ is very important; it allows us to 
design a hypercomplex Fourier transform Tk by prescrib- 
ing the eigenvalues on the basis {'4'j,k,i\ via 

for any set of numbers Xk,fJ.k G C. Indeed, it suffices to 
solve the system of equations 

Xk - \. ( T^ak - sin a ,, /3fc) z^^e'^^'^ 



X / r"+'=~Vo(r-) {^-\Jk+x{'^ e4(^°*")'- dr 
Jo 



r(A + 1) \X + k 
X 



and 



{f{r)xMk) (y) 
X 



X + k + 1 

r(A + i)- 



ak+i + sma 



2{X + k)' 
k + l + 2X 
2{X + k + l) 



k+l 



■k+1 —iak 

I e 



+ 1 + 2A 



ATfcTT"'=+^+'^""2(fc + l + A) 

/' + 00 



with z = (r|y|)/sina, A = (m — 2)/2 and 

2 



to determine the integral kernel K{x, y) in terms of the 
coefficients ak and f3k. 

3.1.3 Inverse transform 

In order to construct the inverse of the general transform 
on the basis {ipj^k.e} we consider the following integral 
transform: 



r(f) (l-e-2»)™/2- 

Proof. The proof goes along similar lines as the proof of 
Theorem 6.4 in |21]- □ 

We then have the following theorem. 

Theorem 1. One has, putting (3q — 0, 



J^Kilp2j,kj) 

2-^ ( X 



if) (y) = P-a,rn / K*ix, y) f (x) dx, 

where the kernel is given by 

K*ix,y) = {A*iw,I) + ixAy)B*{w,z)) 
X e~^^ 

with 

+ 00 



r(A + i) \x + k 



-ak — sm a 



2{X + k) 



■k-ia(k+2j) 



^2j,k.e 



and 



J'Ki^2j+l,kj) 

2-^ f X 



fc + l + 2A „ 
ttfc+i + sma^TT—— ^— -;^/3fc+i 



r(A + l) VA + fc+l 2(A + fc + l 

^, -fc+l -ia(fc+2i + l) / 



A*{w,^=Y.i-\f^k iz)-\Jk+xi^C^iw) 

+ 00 

i?* («;,?) =5](-l)'=+i 6k {l^)-^-\Jk+x{z)Cl+l{w) 
fe=i 

and 7fc,4 e C, z = (|x||y|)/sina, w = (a;, y)/(|x| |y|), 
A= (m-2)/2. 

Similarly as for the transform Fk, we can consecutively 
calculate the radial behavior of the transform J-^ and de- 
termine its action on the basis {^i,k,i}- This yields the 
following result. 
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Theorem 2. One has, putting Sq = 0, 



2-^ ( A 



r(A + 1) VA + /C 



7fe + sm a 



2{X + k) 



Using Theorem [T] and [21 this condition is equivalent with 
the system of equations (fc = 0, 1, . . .) 



A + fc " 2(A + fc) 

7fc + sm a 



X I -A^T-j. + sina— = (-l)^(r(A + l))222^ 



A + /c 



2(A + fc) 



and 



and 



2-^ ( A 



fc + 1 + 2A 



- r(A + 1) V A + + l^'^+i " "2(A + fc + 1) 
Proof. Similar to the proof of Theorem [1] 



4- 



+1 



A fc + 2A ^ 



A + fc 



2(A + fc)' 



X I -A^^. - sin a . ^+ ^^, 4 ) = (-l)'(r(A + 1)^2'^ 



A + fc 



2(A + fc) 



□ 



Combining Theorem [T] and [2l we are now able to con- 
struct the inverse of the general transform J-k on the basis 

Theorem 3. The inverse of J-k on the basis {i/'j.fe.f} is 
given by 

^K^ if) iv) = P-a,m / K{x,y) f{x) dx, 



with 



K{x, y) = {a{w, F) + (x a y)B[w, F)) 



X e 2V 



where 



+ 00 



A{w.^--^^{ak+l3k sina) (z) ^ Jk+x{z)C^{w) 

k=0 



+ 00 



B{w,T) 



and 



Nt = 

k r,2A 



k=l k 



1 



22A(r(A + l))2 VA + fc 

X 



2(A + fc)' 



A . fc + 2A ^ 



A + fc 



2(A + fc)' 



Solving these two equations for 7^ and 5k then yields the 
statement of the theorem. □ 

3.2 Four types of generalized convolution 

3.2.1 Definitions based on an idea of Mustard 

The definitions of the first two types of convolutions are 
based on the observation that in the classical case the 
following interaction between the convolution and the 
Fourier transform holds: 

.F(/*.g) = (2^r/2^(/) J- (5). 

In our case, this leads to the following: 

Definition 1. For f,g G tS(R™) CSCZo.m, the generalized 
convolution f *c.l 9 is defined for x G M™ by 

if *c,L g) {x) := p-^^TJ,^ {Fk (/) (5)) {x). 

Similarly, the generalized convolution f *c,R 9 takes the 
form 

(/ *c,R 9) (x) ■■= p-%,:F^^ [Tk {9) Tk (/)) (x). 

Note that, as already mentioned in the introduction, 

f *c,Rg = 9 *C,L f- 

Taking into account the integral expression for J-k and 
its inverse J^k^, see formula ([8]) and Theorem[3l we obtain 
the following explicit formulas. 

Proposition 2. The generalized convolutions f*c,Lg and 
f *C,R. 9 take the following explicit form: 



Proof. Put K{x,y) = {a{w,z) + {xA y)B{w,z) 
xe~5(™t")(l^l^+l!'l^) where 



if *C,L g) (x) = Ca.rn / / K{u,x) 

JR™ JR"^ JR™ 

X K{t,u)f{t)K{y,u)9{y)dtdydu 



and 



+00 



A{w,7)=Y,{-lf^k iI)~^Jk+xiI)Cj:{w) 

+ CXD 

=^(-l)^-+i4 iT)-^-'Jk+xi^C^lliw) 



if *CM 9) ix) = Ca,m / / K{u,x) 

JR™ JR™ JR™ 

X K{y,u)g{y)K{t,u)f{t)dtdudy, 



with 



k=l 



and with 7^, 4 £ C. We need to have that 

:FK\TKif))^TK{Tj,' (/))-/■ 



Ca,m — Pa,rnP~a,m 



(7r(l-e2-))-^"/^(7r(l 
(27r)-"|sina|-'". 



2mN\-'"/2 
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3.2.2 Definitions using the generalized transla- 
tion operator 

We first define a generalized translation operator related 
to the integral transform Tk defined in Section [3.11 

Definition 2. Let f ^ 5(R"') ® C/o,m- For y G M" the 
generalized translation operator f i — > f is defined by 

^K{T^f){x) = K{y,x) Tk if) (x), X e R". 

It can be expressed, by the inverse of J-k (see Section 
13.1.31) , as an integral operator 

r^fix) = P-c.rn f mTx) K{y, Tk (/) (0 (10) 

Using this generalized translation, we can again define 
two types of convolution for functions with values in the 
Clifford algebra. 

Definition 3. For f,g G 5(IR'") ®Clo,m, the generalized 
convolution f *l g, resp. f g, is defined for x G M™ hy 

if *L g) (x) / [rl'fix)] g{y) dy 

resp. 

{f^H9)ix):^ [ f{y)[T^g{x)\dy. 

Using the integral expression for the generalized trans- 
lation operator (see ([T0|) ): 

rf /(x) = Ca,™ / / K(^x)K{y,u)K{t,u)f{t)dtdu, 

we obtain the explicit formulas for the generalized convo- 
lutions introduced above. 

Proposition 3. The convolutions f *l 9 o,nd f *r g take 
the following explicit form: 

(/ *L g) [x) = Ca.rn i i i K(u, x) 
JR™ jRm JR™ 

X K{y,u)K{t,u)f{t)g{y)dtdudy 

and 

if *R9)ix) ^ Ca^ni / / /(t)if (u, x) 

JR™ JR™ JR"" 

X Kit,u)Kiy,u)giy)dydudt. 

3.2.3 Connection between the four types of con- 
volution and further properties 

For scalar functions the four types of convolution defined 
in the previous subsections are strongly related. 

Proposition 4. Let f,g ^ iS(M'") be scalar functions. 
Then one has 

ft) if *c,L9){x) = if *R9)ix) 



(a) if *c,R 9) ix) = (/ *L 5) [x) 

(Hi) if *L 9) {x) = (g *R f ) ix) 
(iv) if *c\L 9) ix) = (g *c,R f) ix). 

Proof. These relations follow immediately from the ex- 
plicit formulas of the convolutions (see Proposition [2] and 
Proposition |3]). □ 

The following proposition quantifies the difference be- 
tween the left and right versions of the generalized convo- 
lutions. 

Proposition 5. Let f,g& iS(M'") be scalar functions. 
Then one has 

if *L9) ix) - if *R9) ix) 

= if *C,R 9) ix) - if *C.L 9) ix) 

= Ca,m / / / Kiu,x)[Kiy,u),Kit,u)] 

JR'" JR™ JR™ 

X fit)9iy)dtdudy, 

where the commutator of Kiy,u) and Kit,u) takes the 
form 

[Kiy, u),Kit, u)] ^ Biwi,zi)Biw2,Z2) 

^ cota(|i/p + |ip+2|«p) 

X ((y A M)(i Au) - itAu)iyAu)) 

with z\ = (lyllwl)/ sina, zi = (|t||u|)/ sina, w\ — (^1,771), 
W2 = (^j'yi) where y = \y\^i, u = \u\T]iand t = \t\£,2- 

From the definition of the convolutions / *c,l 9 and 
/ *c,R 9 we immediately obtain the following result. 

Theorem 4. For f,g G 5(M™) (g)C?o,m, one has 

^Kif *C,L 9) Palra^K if) (s) 

and 

J'Kif *CM 9) = P»]m^K (5) if) ■ 

Taking into account Proposition 2] we obtain in case of 
scalar functions a similar result for the convolutions f *l9 
and / *R, g. 

Proposition 6. For f,g G 5(M™) scalar functions, one 
has 

^Kif *L 9) = Pa]m^K is) if) 

and 

^liif *R 9) = Pala^K if) (.9) ■ 

3.3 Generalized translation of radial func- 
tions 

The generalized translation defined in formula (jlOp no 
longer equals geometric translation. However, in the spe- 
cial case when we compute the translation of a radial func- 
tion and when a = 7r/2, we still find that generalized 
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translation equals geometric translation. That is the main 
result we will obtain in this section. As a consequence, we 
will be able to give more results on the new convolution 
products when one of the functions is radial. 

The key ingredient we need is a compact formula for the 
integral over the unit sphere 



K{rri,x) K{y,rrj) duj{r]) 



where J^m-i duj{r/) ~ 1. This is derived using the series 
representation of the kernel function and some lemmas 
which are proven in [21]. We state them again for the 
convenience of the reader. 

Lemma 1. Fork,£e N, A = (m-2)/2 andx',y' G S™-\ 
one has 



-1 ~ 



Lemma 2. Fork,£e N, X= (m-2)/2 andx',y' £ S""'^ , 
one has 



(77 Ax') C^^'i(v,x')) (y'Ar;)C^i((77,2/')) 

L 

(fc + l)(fc + l + 2A) ^ , 



with 



+00 



+00 

Ga(w,20 = sina^/3fc iI)-\Jk+x(^C^-lH 
fc=i 

+00 

+00 



k=l 



4A(fc + A + l) 
-4,. ^^^ fe'AyO C^,+\{x\y')). 



We can now show the following result. 



Theorem 5. For x,y £ M™ and r G R+, one has 



K{rrj, x) K{y, rrf) duj{ri) 



= 2^r(A + l)u-^Jx{u) e-*(™t")(l=^l'-|yl') 



withu^ ^^/\x\^ + - 2{x,y) = ^\x-y\ and X 
(to-2)/2. 



Proof. First we slightly rewrite K{x,y) (see Section [3TT 
and K{x,y) (see Theorem[3]) as 

K{x,y)= {Fxiw,I) + (x' A j/)Gx (w, T)) 

K{^y) = {Sx{w,z) + {x: Ay;)Tx{w,^) 

X e^^* 



where Sfc = (/3fc sin a + ak)/N^, tk = -Pk/N^ 
1 



22^(r(A + i))2 \\ + k 



A /c 

ttfe - sma 7:7-^— rrPfe 



2(A + /c)' 



A . fc + 2A ^ 

X I ^ rQfc + Sma ^^^ , n PA; 



A + A: 



2(A + A:)' 



and = x/\x\,y' = w = {x',y'),z = (la;||?/l)/ sina, 

A = (m-2)/2. 

Using these decompositions, we obtain 



K{rri,x) K{y,rri) duj{ri) 



where we calculate the 4 pieces /i, /2, -^3 and I4 separately. 

For Ii , we use the reproducing property of the spherical 
harmonics to obtain 

h = Sx{wi,zi) Fx{w2,Z2) duj{ri) 

00 ^ 

" X! ■^;^^Skak{ziZ2y^Jk+x{zi)Jk+x{z2)C^{{x',y')), 



k=0 



where we use the notations zi — (r|a;|)/ sin a, Z2 — 
(r|y|)/sina, Wi — (77, x') and 1112 — {y',r])- For ^2, we 
use Lemma [T] yielding 



h = I {t] Ax') Tx{wi,zi) Fx{w2,Z2) dLu{T]) 
00 ^ 

sina (x' Ay') ^ ^ ^ tkak{ziZ2T^ 

k=l 

X+l/ii „/\ 



X Jk+x{zi)Jk+x{z2)Cil[{{x\ y')), 
and similarly for 1^ 



h= Sx{wi,zi) {y' At]) Gx{w2,Z2) du!{r]) 

00 ^ 

= - sina {x' Ay') . , , Skl3k{ziZ2)^^ 
— ^-^ A + A; 

k=\ 



X Jk+xizi)Jk+xiz2)C^l,i{x', y')). 
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Finally, we can calculate the term I4 using Lemma [5] as 
follows 

(ryAx') Tx(wi,zi) (y' A 77) Gxiw2,Z2) 



= (sina)2^. 



-tkPk{ziZ2y 



^ Jk+x{zi)Jk+x{z2)Ct{{x\y')) 



fe=i 



X Jfc+A(2:i)Jfc+A(^2)Cfc+]^((a;',y'))- 
Adding these 4 terms then gives 

Ii+h+h+ h 



E 



fc=0 



Sfettfe + (sma) , .. ^fcPfc (2:12:2) 



A + fc 



4A(fc + A) 



-A 



X Jk+x{zi)Jk+x{z2)C^{{x',y')) 



A 



(x' A yOsinaV" . , , fa^fc + Sfc;3fc + sin a tkPk) 
— A + K 



fc=i 



X (zi Z2 ) - Vfe+ A (21 ) Jfe+ A (22 ) q^^i^ ( (x' , y ') ) . 

It is easy to check for all k that tkak+SkPk+sina tkPk = 0, 
so the term in (a;' Ay') vanishes. Similarly we can compute 
that 



A 



A + fc 



. . ,2fc(fc + 2A) 
Skak + {sn.a) ^^^^^ 



tkf^k 



= 2^^(r(A + l))\(A + fc), 



hence we conclude that 

h+l2 + l3 + h 

^ 00 

= 22^(r(A + l))2-^(A + fc)(?iZ2) 



k=0 



X Jk+x{zi)Jk+x{z2)C^{W,y')). 

Now we invoke the addition formula for Bessel functions 
(see e.g. [36], p. 107) yielding 



^JA(^/)-2^r(A)^(A + fc)(ziZ2) 



-A 



fc=0 



X Jk+\{zi)Jk+x{z2)C^{{x',y')) 



with u — ^i^-\/|a;P + lyP — 2{x,y). This completes the 
proof of the theorem. □ 

We can now prove the main theorem in this section. 

Theorem 6. Let f G 5(]R™) be a real-valued radial func- 
tion on M™, i.e. f{x) — fo{\x\) with fo : IR+ 1 — > R, then 

Ty fix) = ^•^i'^yrn/2^-±(cota)(\x\^-\y\^) 



r(A + i)' 

X Hx [Tc (/)] 



\x - y\ 



with A = (to — 2)/2, J-a the fractional version of the clas- 
sical Fourier transform given by the integral transform 



J^a if) iy) = Pa,r. 



^ cot a(|a;| +|y| ) 



fix)dx 



and Hx the Hankel transform defined by 



Hxfis) 



fir) 



Jxirs) 



,2A+1 



dr. 



10 irs)^ 

Proof. If f{x) = /o(|a;|) is real-valued and radial, then by 
means of Proposition [1] we obtain that J^k if) is a radial 
function as well and it coincides (up to a factor) with the 
fractional version J-q of the ordinary Fourier transform. 
More precisely we have 

^Kif) = 

Hence, by definition of we obtain 



-^K(rf /)(x) = 



r(A + i) 



K{y,x) T^{f) (x) 



Taking the inverse and using polar coordinates x = rrj, 
r = we obtain 



fix') = 



(7r(l - e2*"))-"/2 / 



r(A + i) 

X K{y,x) Fa if) i\x\) dx 



r(l 



+00 



-^a (/) (r) 



(r(A + i))2 

X K{rT],x') K{y,r7])duj{T])^ r"^^^dr. 

In view of Theorem [5] this becomes 



ffix') = 



2ao 



(1 



2^Q^)-m/2g-^(cotQ)(|a:'|^-|^/|^) 



r(A + i) 

with u — -^^\x' — v\. Taking into account the definition 

sill Q I ^1 ^ 

of the Hankel transform Hx, we finally obtain 

Ty fix') ^ (1 - g2»a)-m/2g-^(cotQ)(|:cT-|al^) 

\x'-y\ 



X Hx [Ta if)] 



sma 



□ 



In the special case when a — 7r/2, it follows from The- 
orem [6] that the generalized translation operator coin- 
cides, up to a constant, with geometric translation if / is 
a radial function. 

Corollary 1. Let f E 5(M™) be a real-valued radial func- 
tion on R", i.e. f{x) = fai\x\) with /o : IR+ 1 — > M, then 
in case of a = tt/2, one has 

ry^.fix) = ,^^fi\x-y\) 



r(A + i) 



with A ^ (to - 2)/2. 
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Proof. In case of a = 7r/2, the fractional Fourier transform 
Ta reduces to the classical Fourier transform T. Moreover, 
taking into account that for radial functions the classical 
Fourier transform coincides with the Hankel transform Hx 
and that the inverse Hankel transform is given by 



+ 00 



Hxfir) 



Mrs) 



which holds under mild conditions on /, the desired result 
follows. □ 

These results allow us to give more details about the new 
convolution products, when one of the functions involved 
is radial. This is summarized in the following proposition. 



Proposition 7. If g e 5(R™) «) CZo,m and f e is 
a real-valued radial function, then 

if *L 9) = [g *L /) 

= if *R 9) = J^K {9 *R f) 
= Pa!rn^K if) (s) • 

In particular, under these assumptions one has 

f *L 9 ^ 9 "^L f ^ f 9 ^ 9 /• 

Proof. If / is a real-valued radial function, then from 
Proposition[T]we obtain that Fk (/) is a real- valued radial 
function as well. This implies that !Fk if) h = hJ-x if) 
for any Clifford algebra-valued function h. Moreover, from 
Theorem [6] we observe that also f commutes with any 
Clifford algebra- valued function. 

Let us now, for example, show that ^a'(/*l5) = 
Pa\n-^K if) J'K (ff)- The othcr statements are proved sim- 
ilarly. By definition of the transform J^k, the convolution 
*L and the Fubini theorem, we obtain consecutively 

J^K if *L 9) ix) 

= Pa,m I Kit,x) 



[ryfit)]9iy)dyj dt 
Pa.m^^^(i,2:)[Tf/(<)]dtj giy)dy 
Tk ir^f) ix) 9iy)dy 
Kiy,x) Tk if) ix)giy)dy 



= Pa \n^K if) ix) / Kiy,x)9iy)dy 

= Pa]m^K if)ix) Tk ig) ix) , 

where we have used Definition [2] and the fact that 
Kiy,x)TKif)=TKif)Kiy,x). 

Finally, by applying the inverse transform J-j^ on 

J^K if *L 9) = (5 *L f) = Tk if *R 9) = i9 *R f) : 

we obtain 

f "^L 9 ^ 9 "^L f ^ f *r9 ^ 9 *R f- 



4 Convolution products in Ap- 
proach B 

4.1 Definition, eigenvalues and eigenfunc- 
tions 

Let us start by defining the family of transforms we will 
be looking at in this section. 

Definition 4. Denote by I,„ the set {i e Clo^m\i'^ = 
— 1} of geometric square roots of minus one. Let Fi := 
{ii, i^}, F2 := {ifj,+i, ...jim} be two ordered finite sets 
of such square roots, ik G Im,Vfc — l,...,m. The ge- 
ometric Fourier transform (GFT) J-F1.F2 of a function 
f : M™ — )■ Clo^m takes the form: 



•^F,.F.(/)(u) (27r)-^ 



\k=l 



X fix) n 6"*'='='="'= dx. 

\fe=M+l / 

This definition is a special case of the general geometric 
Fourier transforms from [9], where also non- linear func- 
tions in the exponentials are allowed. We will use this 
restriction to guarantee that the inverse transform of any 
GFT is a GFT itself, namely 



-n},{-j™,---,-v+i}- 



Moreover, the restriction also allows us to obtain the eigen- 
values and eigenfunctions of the GFT. They are given in 
the following theorem. 

Theorem 7. The basis {V'ji j2,...j™} 0/ 5(R'") ® CZo,m 
diagonalizcs the GFT. One has 



^Fi.F2ii^juh,...d^) 

= (\{i-^k)Ai^n 



j2,...j™ I n 

k=pi+l 



3k 



Proof. By direct computation, we find 

(27r)'?J^Fi,F2(V'ji,i2,...J„) 
/ M \ 

IM ) 

n ^ ' 



dx 



□ 




dXk 
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using Notation\^ by 
\k=\ ) \fe=p+l / JeJ<?,76{o,i}" 



(1 M 



which is a special case of formula ([6]). 



□ 



ik—l fc— ^+1 k—ni 



with the sign cj ^ - (?wen fey 



An important example of this family of transforms is 
the two-sided quaternionic Fourier transform (qFT). We ^ ^ = (_X)(-J(2*fc+Tfc+i,fc)+i)('5(ii,fc+32,fc+J3,fc>~^\ 

will treat it in detail in Section 1131 ''"'^''^ 



J, 0,7 

k=l 



Remark 4. Note that all subsequent results also hold in where 

the more general setup of functions f : R™ — >■ Clp^q, com- , 
pare However, this is not the case for the transforms g •= < ' ~ 0' 

of approach A, due to the fact that the Dirac operator is \0, if i ^ 0. 

no longer elliptic in arbitrary signature. 

Proof. For the multi-indices j 6 {0,1}'^^'", which we ad 

4.2 Convolution formula based on Mus- '^^^^^ ^ = 1' ^presenting the inte 



tard's idea 



gration variables x, y, z and k = 1, m representing the 
coordinates, we use the notation 



and for = 2, 3 corresponding to y,z analogously. Note 



Based on the idea of Mustard we define a generalized con- 
volution for any geometric Fourier transform. _i^xkUk 1 cos(a;fcUfc), if ji.fe = 0, 

1 - sin(a;fcUfc), if Ji,fc = 1 

Definition 5. For any GFT J-p-^ p^ we define the convo- 
lution *Fi,F2 by 

if *Fr.F2 g){x) {2^T)'^Tp^\p^{TF^Mf)J^F^,F2{9)){x)■ 

Now we want to express the convolution *Fi,F2 by means 
of the standard convolution 

if*g)ix)^ f{y)g{x -y) dy. 



center of the Clifford algebra. We get 



that e^^ fc fe fc jg always real valued and therefore in the 



{2nr{f*F,.P,g){x) 
1 



e 

k—f£ 



(„ /J rn 
k=l fe=/i+l 

/ l[e-'^'^^^g{z) n 



Notation 8. For functions f,g : M'" CZo,m a?*'^ multv 
indices 0, 7 S {0, 1}™ we put 



f^{x) :=/((-l)*^a;i,...,(-l)^'"x„), 



fc=l k=^^+l 



k—m 

The following theorem is our main result. r _i_ JL J!L 

— / TT g«fc2:fo"fc 11 g-ikVkUk j^y^ TT ^-ikVkUk 

Theorem 9. Let J = {0, l}4xm ^i^/i ji,^ + ja.fe + js.fe £ ^=^' ^='^ 

{0, 2} and j^^^k — /or all k = 1, . . . ,m be a set of multi- ^ . . '''^^ 

indices. Any generalized convolution *Fi,F2 from Defini- ^ TT g{z) JJ^ e ^'-'^''^^ JJ^ ^^kXkUk 

tion\^can be expressed as a sum of classical convolutions '^^i *;=/j+i fe=m 
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E 

Je{o,i}3x 

\k=l k=l 



k—fi k—fi j 



"h.k 



n i'^y^'" n 

n ^'^y^'' n ' 
n (*^)''''° n ^jt?"' ^^2^ 



-m k—m 



E 

Jg{o,i}3 
1 



"I /. + 00 



n 



ikXkUk -tkVk^k -ikZkUk 

ii.fc i2,fc j-i.k ^ 



U3,fc 



\fc=l / 

X [n 

\k—m / 



k=fj,+ l 



In order to calculate the integration over the w-variable, 
we use the trigonometric equalities 



sm X sm y sm z 



cos X cos 2/ cos z 



sm X sm y cos z 



- - sm{x + y — z) + sin(j/ + z — x) 
+ sin(z + X ~ y) ~ sm{x + y + z)j 

■ - ^ cos(x + y ~ z) + cos(y + z — x) 
+ cos(z + a; — y) + cos(a; + y + z)j 

: — ^ — cos(x + y — z) + cos{y + z — x) 
+ cos(z + a; — y) — cos(a; + y + z)j 



sin a; cos y cos z = — sin(a; + y — z) — sin(?/ + z — a;) 

+ sin(z + X — y) + sm{x + y + z)^ . 



(11) 



For each index fc = 1, m and ji^k + j2,k + J3,fc = we 



have that the Wfc-integrand takes the form 



^ikXki^k ^—ikUk^k ^—^kZkiJ-k 

jl,k j2.k j3,k 



-<{ cos((-a;fe +yk + Zk)uk) + cos((a;fe ~ yk + Zk)uk) 



+ cos{{xk + yk - Zk)uk) + cos{{xk + yk + Zk)uk) 



and for ji,k + j2,k + h,k = 2 we have 

^-ikVk^k 

-'h.k 



^ikXkiik ^-^kVk^k ^-ik^kii-k 

jl,k j2.k js.k 



cos{{-Xk + yk + zk)uk) 



+ cos((.Tfe -yk + Zk)uk) 

+ (_l)0 3,.+l) ^ _ ) 

- cos((a;fc +yk + Zk)uk) 



Hence for ji_k + J2,fc + J3,fc even we can summarize both in 



^ikXkUk ^-^kVk^k ^-ikZkUk _ 

3l,k 32, k 33, k 



(^ — \y2,k+33.k 



X |(-l)(^'^''=+''^'"i.^+^^>'=+^3..)-l)cos((-.i:fc+yfc + Z,K) 

+ (_l)02..+i)(5a,,,+.,,,+.3,.)~i) cos((a;fc - yk + ZkW) 

^ (_l)03.. + l)(5o,,,+,,,,+,3,,,-l) ^ _ 

+ +^2,^+^3, cos((xfe + yk + Zk)uk) 

with 



1, if£ = 0, 
^0, if £7^0. 

In case of ji_fc + j2,k + j^.k odd, the trigonometric equa- 
tions (HH) immediately show that e*';^^'""'=ej/^'^'="'=ej^'^^'='"= 
equals a sum of sine functions. 
Hence, by splitting the equation 

5{x) f e-""dM 

into real and imaginary part: 



S{x) = — I cos{xu)du, 



= — / sin(a;u) du, 

271" Jr 

one sees that all summands with Ji.fe+j2,fc+j3.fc odd vanish 
after integration over Uk, while the even ones become delta 
distributions. As a result we have 



+00 



JkXkUk ^-tkVkUk -ikZkUk J 



0, if ji,fe + j2,fc + J3,fe is odd, 

lL( — ^ ]j2.k+3a.k V 

X(-l)(^(2*fc+Tfe + l.'^)+l^(''0l,fc+32,fc+J3,fc)"l) 

x,5(a;fc + i-l)*''+^yk + (-iV^+^Zk), 
if ji,k + j2,k + js.fc is even. 
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with j4^fe = Vfc and also 
1 

4'" 



E 



Je{o.i}<'x", 
Vfe=l,...,m:ji fc+j2,fc+j3,fcG{0,2},j4_fc=0 



E 



fiy) 



k=ii+l 



33, k 



\k=l 



X / g{zi,...,Zm) 
m 



fc=i 



X I n (-^fe)''" n 1 '^2/- 
fc=^+i / \k= 



For each index fc = 1, m using 



dix) = / 5(^)'5(2 - x) dz, 



we obtain 

_ 1 
4m 



E E 

Je{o,i}4xm. 0,7G{O,1}" 
Vfc=l,...,m:ji,fc+j2.fe+j3,fce{0,2} j4 fc=0 

X I j^(_l)(-''(2*fc+^fc + l.'') + l)('^Ol,fc+32,fc+j3.fc)"l) 
\fc = l 



1, . . . , m, and obtain 



1 

4m 



E E 

fe=l 



Vfe=l,...,m:ji,fc+J2,fc+J3,fee{0,2}j4,fc=0 

\ 

^-|(j(2*fc+7fc + l,fe) + l)('50i,fc+J2,fc+J3.fc)~^) 

\fe = l 

1 M 



I fc— /i A:— 1 



M+1 



k—m 



□ 



X I n (-*'^)'''' n(-*fc)-"-' 

- (-l)^^yi), . . . , {-lV-{Xra - {-If-Vra)) 

X I n {-ikY'A i\{{iky'A dy. 

\k=fj,+ l J \k=m ) 

Finahy we execute the substitution y'^ — (—1)'* fc = 



Remark 5. One can now introduce an immediate analog 
of the convolution / *2 5 mentioned in the introduction, by 
considering g *Fx,F2 /• 

4.3 Convolution formula based on gener- 
alized translation 

We define a generahzed translation operator related to the 
GFT J-p-^ . Contrary to our original definition in formula 
, we now have to take into account that the kernel con- 
sists of two parts. 

Definition 6. For any GFT J-Pi.F^ general transla- 
tion operator Ty^'^^ is defined by the relation 

FF,Mr^^'^^f){u) 

■■= (f[e'"''''Aj'F,Mmu) f fl ' 

With calculations analogous to the ones in the previous 
section we can express the generalized translation operator 
Ty'^'^^ by means of the standard translation 

'Tyfix) f{x ~ y). 

Theorem 10. Let J = {0, 1}-*^" with ji,fc + j2,fc + j3,fc £ 
{0,2} and j^^k — for all k — l,...,m be a set of 
multi-indices. The generalized translation operator Ty^ '^^ 
from Definition \d\ can be expressed as the sum of classical 
translations Tyff^{x) — f^{x — y*^) using Notation\^ and 

((-l)^^yi,-.-,(-l)"*'"2/m) hy 



c- 7 - 



jeJ cl>ne{o,i}'-' 



1 M A* 

X n(^fe)''-' n iiHfeF'''" 

k—fi k—1 k—1 

X Tyff'^ix) 

m /i+l 



X 



k=n+l 



k=fi+l 
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Proof. Similar to the proof of Theorem [3] 



□ 



Using this result, we can now give an explicit expression 
for the convolution product defined using the translation 
operator. 

Corollary 2. Let J = {0, Ij^^" with ji,fc + j2,k + js^k e 
{0, 2} and j^^k — for all k — 1, ... ,m he a set of multi- 
indices. The convolution defined by 

{f*F,.F.9){x):= [ fiy)[ry^^-^^9{x)]dy 



with Ty'-'^^ from Definition\E\ can be expressed 
of classical convolutions using Notation\^ by 

=—y y c,^-. 

JeJ0,7G{O,l}™ 

f^YliikV'^' HHkY'-' n^-*^)'''' 



as a sum 



k=fi 



fe=l 



k=l 



* 9^ n (-^fe)'''" n (-^fe)'''" n ( 



(x). 



k=fi+l 



k=ii+l 



Note that a similar expression as in Corollary [2] can be 
obtained for 



T-Fl,-F2 

y 



f{x)] g{y) dy. 



4.4 Special case: the two-sided quater- 
nionic Fourier transform 

4.4.1 Definition of the qFT 

The quaternion algebra III is isomorphic with the Clifford 
algebra CZo.2 under the identification i = ei, j =62 and 
k = 6162. 

Let /i, 1/ G H be quaternions with = = —1. Then, 
following 05], we define the two-sided qFT as 

^'''''i.f){yi,y2) 

:=(27r)-i / e-^''■'^y'f{xl,X2)e-'"=^y^dxldx2 



for functions / G Li(R^;H) where we have introduced a 
different normalization (27r)~^. The first definition of this 
two-sided transform, with = j and v = k., was intro- 
duced in the Ph.D. thesis ^7}, see also [52]. In earlier 
work, a one sided version was given by Ernst et al. [30] 
and by Delsuc [23], although these authors use an adap- 
tation of the quaternion algebra. The applicability of the 
qFT to color image processing was first demonstrated in 
[33] using a discrete version. At that point, the switch 
was made to two general orthogonal axes fi and v instead 
of j and k. Indeed, for color image processing there is an 
arbitrary but preferred axis of the grey-line in the color 



space, so the transform kernel axes are generally aligned 
to or perpendicular to this axis. At the same time ( |44| ) . 
a change was again made to one-sided transforms, mostly 
driven by the complexity of the resulting operational for- 
mula when using the two-sided qFT definition. Finally, 
the orthogonality condition on /i and v was relaxed in 
|35| . For a recent review on the use of the qFT in image 
processing, we refer the reader to 

4.4.2 Convolution for the qFT 

In this section, we discuss the Mustard and generalized 
translation definitions of the convolution product for the 
qFT. First, we give the interaction of the qFT with the 
ordinary convolution. To that aim, we need the following 
definition. 

Definition 7. For an invertible multivector b G Clo^m o,nd 
an arbitrary multivector a G ClQ_rn we define the commu- 
tative and anticommutative part of a with respect to b by 

acO{b) =^(a + b'^ab), 

"ci(fc) =^(a - b^^ab). 

Using this definition, the quaternionic Fourier transform 
of the ordinary convolution (/ * g){x) is given by 

T^-''{f^g)iu) 
1 

j,k=0 

This formula is a special case of the convolution theorem 
for a general GFT, obtained in [TU]. In the discrete case, 
a similar formula was earlier obtained in [29J. Note that, 
in particular. 



J-^^^ if*g)^27rT'^'''{f)T^^''{g). 



(12) 



The Mustard convolution, as obtained explicitly in The- 
orem [5] for an arbitrary GFT, reduces in the case of the 
two-sided qFT to 



{f*q9){x) = ^ 



E 



Je{o,i}*x2, 
Vfe=l,2:ji,fc+i2.fc-l-i3,fc6{0,2}j4,fe=0 

* {{-y.y''-'g^[-vY''^{vy'-^) {x). 

It can be checked that this expression coincides with the 
more symmetric formula 

1 ^ ^ 
{f*q9){x)^^ y^ y^ Cn,32MM 

ji,j2=0 ki,k2=a 

X {{fi^'f^'v^') * {n''g'''i^^^)){x) 
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where 



References 



and 

/'Ha;i,a;2) = /(xi,(-l)'=^X2), A:i e {0, 1} 
g^Hxi,X2) = g{{-lf^xi,X2), fc2 e {0,1}. 

The advantage of the Mustard definition is that in this 
case, contrary to formula (|12l) . the fohowing holds: 

The situation is quite different for the convolution de- 
fined using the generalized translation. First of all, it can 
easily be proven (see e.g. [14]) for any function / that the 
generalized translation in the case of the qFT coincides 
with geometric translation: 

As a consequence, the associated convolution also coin- 
cides with the classical convolution, i.e. 

(/ *l^.g){x) = {f*g){x). 

This can also be proven starting from the result in Corol- 
lary H 

5 Conclusions 

In this paper, we have studied two conceptual ways of 
defining convolution products, namely using the method 
of Mustard and using the generalized translation opera- 
tor. We applied these ideas to two important families of 
hypercomplex Fourier transforms. A summary of our re- 
sults can be found in Table [T] 

We expect that in particular the Mustard convolution 
will find many applications in color image processing. Cur- 
rently, the design of filters for such images using hypercom- 
plex methods is hindered by the lack of a proper convolu- 
tion theorem. Our results now enable the development of 
a complete theory of linear system design in the quater- 
nionic case (as well as in higher dimensions, which may 
be interesting for other applications). A Fourier domain 
analysis, using the Mustard convolution, of the color edge 
filter constructed in [32] will serve as the guiding example 
to achieve this. 
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Table 1: Results on convolutions for the two approaches 
A and B 





Approach A 
eigenfunction approach 


Approach B 
generalized roots of —1 approach 


Definitfon 

Eigenfunctions 
Eigenvalues 


formulas d?]), % 

complex numbers 
Theorem [T] 


Definition H 

?/; ■ ■ 

elements from Clo.m 
Theorem [7] 


Mustard convolution 
Expression for convolution 


Definition [T] 
Proposition [2] 


Definition [5] 
Theorem [S] 


Generalized translation 

action on: 

Translation convolution 
Expression for convolution 


y 

Definition [2] 
radial functions 
Theorem [6] 
Definition [3] 
Proposition [3] 


tFi.F2 
V 

Definition [6] 
arbitrary functions 
Theorem [To] 
CoroUaryU 
Corollary H 
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